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Abstract 

cn 

^ We show that the dark matter (DM) could be a hght composite scalar rj, emerging 

00 

^ from a TeV-scale strongly-coupled sector as a pseudo Nambu-Goldstone boson (pNGB). 

Such state arises naturally in scenarios where the Higgs is also a composite pNGB, as in 
^ 0(6)/0(5) models, which are particularly predictive, since the low-energy interactions 

^ of rj are determined by symmetry considerations. We identify the region of parameters 

CN where 77 has the required DM relic density, satisfying at the same time the constraints 

^T^ from Higgs searches at the LHC, as well as DM direct searches. Compositeness, in 

addition to justify the lightness of the scalars, can enhance the DM scattering rates 
^ and lead to an excellent discovery prospect for the near future. For a Higgs mass 

d mh — 125 GeV and a pNGB characteristic scale / < 1 TeV, we find that the DM mass 

is cither ~ 50 — 70 GeV, with DM annihilations driven by the Higgs resonance, 
or in the range 100 — 500 GeV, where the DM derivative interaction with the Higgs 
becomes dominant. In the former case the invisible Higgs decay to two DM particles 
could weaken the LHC Higgs signal. 



1 



1 Motivation 



From a theoretical point of view light scalar particles are unnatural, unless a suitable structure 
protects their mass from large quantum corrections. In the Standard Model (SM) there is a 
compelling case for a light scalar, the Higgs boson, and many efforts have been made to address 
the associated hierarchy problem. The dark matter (DM) energy density of the Universe could 
also be accounted for by a new light scalar. The minimal realization, extensively studied in 
the literature (see e.g. Ref. [1]), consists in adding to the SM a gauge singlet real scalar 77, 
and assuming that it is stable due to a parity rj — )■ —17. The model is quite predictive since it 
only depends on two extra parameters: the singlet mass and its "portal" coupling to the 
Higgs boson, A. This apparent simplicity, however, calls for an ultraviolet completion. As we 
said, light scalars are unnatural in quantum field theories, unless they are accompanied by 
new ingredients, such as for example supersymmetry, and these deeply affect the dynamics of 
the DM. 

Light scalars can also be natural if they are not elementary particles, rather composite 
states emerging from a strongly-coupled fundamental theory. Similarly as pions in QCD, light 
scalars can appear as pseudo Nambu-Goldstone bosons (pNGBs) arising from the spontaneous 
breaking of the global symmetries of the strong sector around the TeV scale. A well-known 
example is the symmetry breaking pattern SO{5)/SO{4), that generates four pNGBs with 
the quantum numbers of the SM Higgs doublet H, the weak gauge group being embedded in 
5*0(4) ~ SU{2)l X SU{2)ji [2]. The most minimal extension to this symmetry pattern is the 
S0{6)/ S0{5) model, that contains five pNGBs in the spectrum: the Higgs doublet H and a 
gauge singlet r] [3]. Interestingly, the pattern S0(6)/ S0{5) is the minimal example with a 
known ultraviolet completion in terms of techni-quarks [4] . 

The purpose of this article is to study under which conditions rj is a suitable DM candidate. 
The properties of 77, being a pNGB, are determined by the global symmetries of the strong 
sector, and by the way these symmetries are explicitly broken by the couplings of the SM 
fields to the strong sector. We will analyze the phenomenology of this composite DM particle 
and its interplay with the composite Higgs doublet H. 

The DM couplings are substantially different from the non-composite singlet case. This is 
mainly due to new non-renormalizable interactions between 77 and the SM fields, arising from 
operators of dimension-six suppressed by 1/p where / ~ TeV is the decay constant of the 
NGBs. There are two types of these interactions: (i) derivative couplings between 77 and H, 
fully determined by the S0{6)/ S0{5) structure, that scale as / P where p is the relevant 
momentum in a given process, and (ii) direct couplings between r] and the SM fermions, 
arising from the explicit breaking of the global 5*0(6) symmetry, that scale as m^p/ p where 

is the mass of the fermion. These interactions make the phenomenology of the composite 
scalar DM substantially different from that of an elementary scalar (similar candidates of 
composite scalar DM have been also proposed, for instance, in the context of techni-colour 
theories [5] and gauge-Higgs unified models [6]; see also Ref. [7]). 

In section 2 we present the effective lagrangian for the composite 7], defined from symmetry 
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considerations and naive dimensional analysis (NDA). In section 3 we compute the composite 
DM relic density. In sections 4 and 5 we study respectively the constraints coming from the 
Higgs searches at the LHC and the DM direct detection experiments. We combine these results 
in section 6, where we identify the regions of parameters that fulfill all the phenomenological 
requirements. We conclude in section 7 and leave the technical details of the composite models 
for appendix A, and the lengthy expressions for the relic density and direct detection cross 
sections for appendices B and C, respectively. 

2 Light scalar dark matter from a composite sector 

We consider theories with a light scalar sector consisting only of a Higgs doublet H and a 
singlet T) with parity 



that makes r] stable. The lightness of these states is a consequence of a global symmetry. The 
simplest realization consists of having a new strong sector with a global symmetry breaking 
pattern given by 0(6) — t- 0(5). Five NGBs emerge from this spontaneous breaking: a 
5 of 5*0(5) which decomposes as 4 © 1 ~ (2,2) © (1,1) under the subgroup 5*0(4) ~ 
SU (2) L X SU {2)r. The 4 is identified with if, while the gauge singlet is rj (see appendix A for 
details). Thinking of 50(6) as the rotation group in a six-dimensional space, the symmetry 
under which r] shifts corresponds to a rotation in the 5-6 plane, 5'0(2)5„e = 5*0(2).^ ^ U{1)^, 
while Eq. (2.1) corresponds to the six-dimensional parity of 0(6), = diag(l, 1, 1, 1, —1, 1). 

If it were for the strong sector alone, the pNGBs H and r] would be massless and would 
only interact derivatively. However, the SM gauge bosons and fermions couple to the strong 
sector breaking explicitly 0(6), generating non-derivative interactions between rj and the 
SM fields. We will assume that is preserved by the SM couplings to the strong sector. At 
the one-loop level these couplings induce a potential for H and rj, that is eventually responsible 
for electroweak symmetry breaking (EWSB), with v = \/2{H) = 246 GeV and {r]) = 0. 

This scenario, apart from giving a solution to the hierarchy problem, since the Higgs is 
naturally a light state, also provides a good candidate for DM, the extra state r], that is also 
naturally light. At energies below the strong scale, denoted by mp ~ few TeV, the lagrangian 
for rj at the lowest order in a 'rf / P expansion, where / is the pNGB decay constant, is given 



r] ^ -r] , 



(2.1) 



by 




(2.2) 



+ 




(2.3) 



where 



V{r^,H) = ^fiy + \\H\'r^' + -- - . 



(2.4) 
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For simplicity, only the interactions with the third family of SM quarks are shown, that are 
the ones that play the most prominent role. Let us remark that these operators are not 
suppressed by the mass of the heavy composite states rrip but rather by the smaller scale 
/ ~ rup/gp where gp is the inter-composite coupling, expected to satisfy 1 ^ (^p ^ 47r. 

A few comments on the lagrangian in Eqs. (2.2)-(2.4) are in order. The only interaction 
that preserves the Nambu-Goldstone shift symmetries is the derivative term in Eq. (2.2) that, 
as shown in appendix A, is fully determined by the S0{6)/ S0{5) symmetry. The coefficients 
Ct^b, instead, depend crucially on how the fermions couple to the strong sector and break 
explicitly the SO {6) symmetry. In general, we expect q (cb) to be an 0(1) complex number 
whenever the U{l)rj symmetry is broken by the and/or by the tji (6/j) coupling to the 
strong sector; otherwise they are zero. Both possibilities are realized in simple models, as 
described in appendix A. Notice that the interactions of Eq. (2.3) can also be induced from 
operators of the form 

qlYdpgL"^ , (2.5) 

and similarly for tji and fo/j. Indeed, using the equations of motion for the fermion we can 
rewrite the above operator as those in Eq. (2.3). 

Let us briefly discuss the potential V{ri,H). It can only arise from loop effects involving 
the SM fields. Since the SM gauge interactions preserve the U{l)n symmetry, /^^ and A can 
only be generated from fermion interactions that break U{l)j^. Since A is further protected 
by the symmetry under which H shifts, we generically expect 

A < fil/f ■ (2.6) 

Then, the mass of t] is given by = /i^ + Af ^ ~ /i^, since electroweak precision measurements 
require [8] 

It is easy to construct models where the U{l)r, symmetry is either broken by the top or the 
bottom coupling to the strong sector, giving a one-loop mass for t] of the order of 500 or 
50 GeV respectively, as shown in appendix A. We also notice that in the potential we have 
ignored terms beyond the quadratic order in the t] field, since they do not play any important 
role in our DM analysis. 

Apart from Eq. (2.2) and Eq. (2.3), there is also the possibility to have at order rf / P the 
interaction terms ^ 

^ [cF,Fp,F^- + CF.Fp^F^''^ . (2.8) 

•' F=B,W,G 

These couplings do not respect the U{l)r, symmetry and therefore can only be induced by 
loops involving heavy composite states that see the f/(l)^-breaking through mixing with the 
SM states. We estimate them to be Ci?^, Cpr) ~ {Qf I QpYij^r^/'^p)'^ i smaller than contributions 
coming from SM loops. We will neglect them from now on. 
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In summary, the interactions between our DM candidate rj and the SM particles are 
controlled by four parameters: Cf,, Ct, f and A. For our DM analysis below we shall consider 
the following values for these parameters. Since the bottom quark plays a main role in the 
computation of the DM relic density, in order to highlight the effects of compositeness we will 
assume that its couplings to the strong sector break f/(l)r?, so that Cf, 7^ 0. More specifically, 
motivated by the models described in appendix A, we will consider the range 



where small (large) values of a and b correspond to small (large) couplings of bi^R to the 
strong sector. For the top coupling q we will consider two cases: 

• Case 1: U{l)n is broken by the top and we fix for definiteness q = 1/2 as in minimal 
models (see appendix A). The top loops make 77 generically heavier than the Higgs. 

• Case 2: The top interactions do not break f/(l),,. We then have Q = and does 
not receive contributions from top loops. In this scenario the DM particle rj can be 
naturally lighter than the Higgs. 

For the decay constant / we will take / = 500 GeV and / = 1 TeV. Recall that although 
composite Higgs models generically predict / ~ f ~ 246 GeV, electroweak precision mea- 
surements require Eq. (2.7) to hold. Finally, we will leave A and as free parameters to be 
determined by the requirement of 77 to be a realistic DM candidate. 

The couplings of t] to the 1st and 2nd family quarks will be important only when con- 
sidering direct DM detection. If not otherwise specified, we will assume for simplicity family 
independent coefficients, = Cc = Q and Q = = c^. This choice is partly motivated by 
the bounds on flavour violation [3]. 

To recover the predictions of the renormalizable model [1] , that we will refer to as the 'non- 
composite case', we can take the limit / — )• 00. In this case only A controls the interactions 
of r] with the SM fields. 

3 Relic density of the composite dark matter 

In the standard cosmological framework, DM was kept in thermal and chemical equilibrium 
as a consequence of its interactions with the SM particles. As the Universe expanded and 
cooled, the number density of DM particles decreased until they could not annihilate any- 
more, freezing out by the primordial thermodynamical equilibrium. Thenceforth their number 
density remained constant, fixed to the value observed today. This picture is described by a 
Boltzmann equation that, under certain assumptions, has the useful approximated solution [9] 



+ a + ib , 



a, be [0, 1] , 



(2.9) 



3 ■ lO-^^cm^s-i 



(3.1) 
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where = p^/ pc is the ratio between the energy density of DM and the critical energy 
density of the Universe, h = Hq/ (100 km s^^Mpc^^) is the reduced value of the present Hubble 
parameter, and (crvrei) is the thermal-average of the total annihilation cross section of DM 
particles, times their relative velocity. The present experimental value is VLoMh? = 0.1126 ± 
0.0036 [10]} 

Although the singlet pNGB rj interacts strongly with the heavy resonances of the composite 
sector, its interactions with the SM, including the composite Higgs, are weak, as shown in 
Eqs. (2.2)-(2.4); as a consequence, its annihilation cross section into fermions, gauge bosons 
and the Higgs falls into the ballpark suggested by Eq. (3.1). In appendix B we display all 
the relevant annihilation cross-sections that enter in the Boltzmann equation, both in the 
non-composite and in the composite case. The Feynman rules needed in the computation are 
collected in Table 2. 

The numerical solution of the Boltzmann equation is shown in Fig. 1 as a function of 
rriri, for two representative values of A. To understand the behavior of the relic density with 
rriri, one should keep in mind that each annihilation channel rjr] <H- XX opens only for t] 
masses above the threshold value mx\/l — v^^jA, with Vrei ~ 2/3 at freeze-out. In order to 
analyze the effect of the composite interactions, let us recall that in the non-composite case 
the annihilation cross section is entirely determined by the value of A (see Eqs. (B.13-B.16)). 
In the composite case the situation is substantially different. One can identify four interesting 
mass regions: 

1. The low-mass region, < 50GeV. The annihilation cross-section is dominated by the 
hh channel, that is enhanced by the direct DM— bottom coupling proportional to cj,, 
shown in Eq. (2.3). As a consequence, one can reproduce the DM relic density even for 
values of A smaller than in the non-composite case. 

2. The resonant region, ~ vrthj^. The presence of the Higgs resonance enhances the 
annihilation cross section by several orders of magnitude. In the non-composite case 
this is the only possibility to fit the DM relic density when A is relatively small, however 
for too small values (A < 10"'' at rrih = 125 GeV) even this enhancement is not sufficient 
and the model is excluded. On the contrary, in the composite case one can reproduce 
the DM relic density even for vanishing A, because of the extra derivative contribution 
to the rj-rj-h vertex, coming from Eq. (2.2). 

3. The cancellation region, ~ A/^/2. A cancellation can occur between the derivative 
and the A contributions to the rj-rj-h vertex (see Eq. (2.2) and Eq. (2.4), respectively), 

^ An alternative mechanism to generate the DM rehc density is freeze-in, that is reahzed when the DM 
interactions with the SM are so weak that the DM species was never in equihbrium. In this case SM particles 
can slowly annihilate (or decay) into DM generating the required relic density. For the case of a scalar singlet 
T] coupled to the Higgs portal, freeze-in is possible both for of the order of the electroweak scale [11] and 
for rrijj < 1 GeV [12], as long as A < 10"^*^. Here we will consider only composite models that induce a 
larger value for A as well as a large DM-Higgs derivative interaction. Therefore we will focus on the freeze-out 
scenario. 
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A = 10-', = 125 GeV A = 10'^ m,, = 125 GeV 




50 100 150 200 50 100 150 200 

m,, [GeV] [GeV] 



Figure 1: Relic density of the scalar singlet DM as a function of its mass m^^. We take 
ruh = 125 GeV and X = 10^^ AO '^y' ^^e left (right) panel. The dashed curves correspond to 
the non-composite case, the red (blue) curves to the composite case with / = 1 TeV (500 GeV). 
The red and blue bands describe the variation of the DM— bottom coupling Cf, in Eq. (2.9), 
between a = b = (dark lines) and a = b = 1 (light lines). Above the threshold for annihilation 
into ti, the result depends also on the value of the DM— top coupling: the dotted and solid 
lines correspond to Ct = 1/2 (Case 1) and q = (Case 2), respectively. 



that may suppress the annihilation with s-channel Higgs-exchange when this is too large. 
The cancellation condition reads 

s = 2\f{l-i) , (3.2) 

where s = 4m^/(l — v'^^Ji). For example, for A = 10~^ and / = 500 GeV, taking 
into account the freeze-out temperature one finds that the effect of the cancellation is 
maximal for 70 GeV < < 100 GeV, as one can see in the left panel of Fig.l. In 
the non-composite case (/ — i- oo), when A is sufficiently large the relic density is too 
suppressed for all values of above the resonant region. In the case of compositeness, 
instead, the cancellation enhances the value of and can make it compatible with the 
DM relic density. For larger /, the cancellation region moves out to larger values of m^, 
where the tt channel opens. If the top has no direct coupling to rj (solid lines in Fig. 1), 
the cancellation remains effective. If instead the top couples directly to DM (dotted 
lines), this reduces the relic density already for values of in the cancellation region. 

4. The high-mass region. As the DM mass increases, ^ and ^ A/^/2, the 

Tj-Tj-h vertex is more and more dominated by the derivative coupling of Eq. (2.2), whose 
strength is uniquely fixed by /. As a consequence, the annihilation rate through the 
Higgs portal becomes larger (typically too large for / = 500 GeV). Still, for / < 1 TeV 
the annihilation rate is of the correct order of magnitude, as can be seen in the right 
panel of Fig. 1. The proper relic density can be reproduced even for a vanishing A. 

The analysis of the rj relic density will be completed in section 6, where we will show the 
contours for Q^] = ^dm in the plane (m^. A). Before that, we need to discuss the constraints 
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coming from LHC Higgs searches (section 4), and from DM direct detection experiments 
(section 5). 

4 Constraints from Higgs searches at the LHC 

The SM Higgs boson mass was restricted to be heavier than 114 GeV by LEP measurements. 
This hmit apphes also to a Higgs boson that decays invisibly, since even in this case the bound 
from associated production with a Z boson holds [13]. The LHC 95% C.L. upper bound in the 
low mass region is presently 127 (129) GeV from the CMS (ATLAS) experiment, assuming 
a SM-like Higgs. Here we will not discuss the very high mass region [mh ^ 600 GeV) that 
is still allowed.^ In the low mass region there is also a hint for a Higgs with m/j ~ 125 GeV, 
coming mostly from the decay channel /i — )■ 77 [14]. 

In the singlet extension of the SM, invisible Higgs decays are allowed as long as mh > 2m^. 
The observation of the Higgs boson would imply an upper bound on the invisible decay width. 
Tiny = r(^ VV)- Vice versa, a signal suppression could be explained by a sizable Tiny. In 
our scenario we find 

In addition to the invisible decay channel, compositeness implies further deviations from the 
SM predictions, because the Higgs couplings to gauge bosons and fermions are modified at or- 
der ^ [8] , as we will describe below. These modifications affect both the Higgs production and 
the relative branching ratios for the Higgs decays into visible channels. Reduced (enhanced) 
couplings would weaken (strengthen) the Higgs signal at the LHC. 

Let us discuss first the non-composite scenario, when ^ vanishes. In this case the LHC 
exclusion limits on the Higgs mass are conveniently described by the ratio 

asAiipp ^ h SM) dsMiPP h) TsM + T^inv 

where Tsm is the total Higgs width in the SM, and Tiny is given by Eq. (4.1) in the hmit 
^ — )■ 0. Note that the Higgs production cross section is unchanged w.r.t. the SM, however 
the Higgs visible branching ratio becomes smaller than one, as long as ruh > 2m^. Actually 
for |A| > 0.1 the channel h ^ rjr] dominates. 

In the case of composite h and r], there are several important differences. First, there is 
another h — i] — t] coupling besides A, that arises from the derivative interaction in Eq. (2.2). 
As a consequence, the two contributions to the decay amplitude in Eq. (4.1) can cancel each 
other for A > 0, or add up for A < 0. In the former case the Higgs signal at the LHC is 
maximal for A ^ m^/2/^, rather than for A = 0. 

^ We have checked that, even for a very heavy Higgs boson, with m/j > 600 GeV, the composite singlet 
can still provide a good DM candidate. 
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Second, as we mentioned, there are order ^ modifications of the Higgs couphngs, that are 
specific to the composite nature of h independently from the existence of a hght singlet rj. 
For concreteness, we will center here on the SO{Q) / SO{b) models described in appendix A, 
assuming Case 2 for the top quark. The relevant Higgs couplings are modified as follows: 



9hvv = 9hvvV^ for V = W,Z, (4.3) 
9hu = ^/fifv^W, (4.4) 

9hbb = 9lb{ yf^T^ ■ ^^'^^ 

The Higgs coupling to gluons, which is crucial for the total production cross section, is mainly 
generated by top loops involving the coupling ghtt linearly. Therefore we have'^ 

9r..^9irM = 9Ul^^- (4-6) 

9htt 

Similarly to the case of gluons, the Higgs coupling to photons is dominated by a top quark 
loop and a W^-loop, that involve linearly ghtt and 9hww respectively. Since both couplings are 
corrected in the same way (see Eqs. (4.3) and (4.4)), one finds 

9h,, ^ 9'h^',Vi^^ ■ (4.7) 

In order to compare these effects of compositeness with the LHC results, one should 
notice that the most sensitive channels in the present analyses by ATLAS and CMS are 
h — )• 77, WW, ZZ. Therefore, what is actually measured is not the total Higgs signal strength 
/i, rather the signal strength in the 'gauge' decay channels only. As a consequence, the relevant 
quantity is 

a{j)p^h^VV) _ a{pp^h) T{h^VV) Tsm 
crsuipp VV) asuipp -> h) Tsuih VV) V^^^ + Tinv 

^ (1 - 0(1 - 0-7 'SM + ^SM ^4_g) 

rS(i- 207(1 -0 + r^]:.(i-0 + r... 

where V = •y, Z, W, and we assumed that the Higgs production is dominated by the gluon 
fusion channel, controlled by Eq. (4.6). Here ^'sm^ is the total visible width in the compos- 
ite scenario, while T'plj and T^m are the total SM widths into fermions and gauge bosons, 
respectively. Note that, in the non-composite case, Eq. (4.8) reduces to Eq. (4.2). 

In Figs. 2,3 we display our prediction for ^vv as a function of A. Since the gauge channels 
are the most sensitive at the LHC, the experimental constraints on /i [14] are more appropri- 
ately interpreted as constraints on fivv- The red shaded regions in Figs. 2,3 are disfavoured 
at 95% C.L.. Let us discuss the 'light' and 'heavy' Higgs scenarios in turn: 



One should remark that, in composite Higgs model, the suppression of the top quark contribution to 
dhgg could be compensated by the contributions of heavier states, such as vector-like fermions that accompany 
the top, or other resonances of the strongly-interacting sector. However, in most cases it was found that no 
enhancement of the Higgs production occurs, but exceptions are possible [15]. 
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Figure 2: The strength fiyv of the Higgs signal in the gauge channels, defined in Eq. (4-8), 
as a function of the h — r] quartic coupling X (negative in the left panel, positive in the right 
panel). We chose rrih = 125 GeV and two values for the scalar singlet mass rrir,, one larger 
than mh/2 and the other smaller. The dotted curves correspond to the non-composite case. 
The dashed (solid) curves correspond to compositeness with / = 1 TeV (500 GeV). We took 
into account the order ^ corrections to the Higgs couplings to vector bosons and fermions, see 
Eqs. (4.3)-(4.7). The shaded region is disfavoured by the LHC Higgs searches (we roughly 
extracted the 95% C.L. lower bound on fiyv from Ref. [14])- 




Figure 3: The same as in Fig. 2, except for rrih = 145 GeV. For this mass ATLAS and GMS 
put an upper bound on fiyv- 



• SM-Uke Higgs with rrih < 130 GeV (Fig. 2) 

If the signal around ruh = 125 GeV were confirmed, the ratio ^vv should lie relatively 
close to one. More precisely, the CMS (ATLAS) collaboration finds the largest excess at 
rUh = 124(126) GeV, with /iyy = 0.6 - 1.2 (0.6 - 1.3) at la [14]. Given this large statistical 
uncertainty, in Fig. 2 we roughly estimated the 95% C.L. lower bound as fi> 0.3. 

In the non-composite case, for > rrih/ 2 one has /iyy = 1, while for m,, < rrih/ 2., in 
order to avoid a too large Tinv, one is forced to take sufficiently small values of A. 

In the composite case instead, when m.^ > mh/2 we find that /lyy can be larger than 
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one. Indeed Eq. (4.8) shows that, even though the Higgs production is suppressed by a fac- 
tor (1 — the gauge decay channels receive a (1 — ^)^/(l — 2^)^ enhancement relatively 
to the fermion channels (which dominate the branching ratio up to rrih ^135 GeV). As a 
consequence, one finds that fiyv can be as large as 1.2 for / as small as 500 GeV. When 
rriri < mh/2, the signal is suppressed by Tinv However larger values of A are allowed w.r.t. 
the non-composite shown in the right panel of Fig. 2, because of the cancellation 

between the two terms in brackets in Eq. (4.1). 

• Invisible Higgs with rrih > 130 GeV (Fig. 3) 

In this region the LHC searches exclude a SM-like Higgs, by requiring the ratio ^vv to 
be smaller than one. In our scenario, this can be accounted for by the combined effect of 
invisible decays and reduced SM couplings, as illustrated in Fig. 3, where we fix m/^ = 145 
GeV (for this value of the Higgs mass the most sensitive channels are h — )■ WW — 2/2z/ and 
h — )■ ZZ — )■ 4/). In this case the 95% C.L. upper bound from ATLAS and CMS is roughly 
f^vv ^ 0.4. There are wide regions of parameter space that survive the LHC constraint, as 
shown by the curves lying in the unshaded region in Fig. 3. 

5 Constraints from dark matter searches 

The direct detection of DM is based on the observation of the elastic scattering between non- 
relativistic DM particles from our galaxy halo and cryogenic nuclei in targets, that produces 
a nuclear recoil. In our scenario the interaction of the composite DM rj with nucleons arises 
through the direct couplings Cq (g = u,d,...) between DM and quarks, given in Eq. (2.3), 
and through the t-channel exchange of a Higgs boson, which in turn couples with quarks; in 
both cases heavy quark loops also induce a coupling to the gluons in the nucleons. These 
processes give rise to a spin-independent elastic cross section asi, potentially within the reach 
of present and future experiments. Note that Im(cg) contributes only to the spin-dependent 
cross-section asn, which is always much smaller than cxsi (the latter being enhanced by the 
coherent interactions on protons and neutrons in the large target nucleus). 

The XENONIOO experiment recently set the best experimental upper limit on asi [16], and 
it plans to considerably improve its sensitivity already within the end of 2012. The excluded 
region at 90 % C.L. is shaded in green in Fig. 4. On the other hand, the DAMA/LIBRA 
collaboration (confirming previous results by DAMA/Nal) has provided the first evidence of 
an annual modulation in the event rate, that could be due to DM, with 8.9 standard deviations 
from the expected background [17]. Similar signals have been observed by CoGeNT [18] 
and by CRESST-II [19], which both measure an excess of low-energy nuclear recoil events. 
CoGeNT also observes a seasonal variation. The analysis performed e.g. in Ref. [20] shows 
that a DM candidate with moM — (10 — 20) GeV and asi ~ (1 — 3) • lO^'*^ cm^ can account 
for the excess events reported by each of these experiments. This favoured region is shaded 
in orange in Fig. 4. One can appreciate the well-known tension between these signals and the 
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A = 10-', in,, = 125 GeV A = 10'^ m,, = 125 GeV 




wi,, [GeV] OT^ [GeV] 

Figure 4: The spin-independent cross section for the elastic scattering of the DM candidate r) 
off nuclei. The green shaded region is excluded by XENONIOO [16], while the orange shaded 
region roughly corresponds to the excess events reported by DAMA, CoGeNT and CRESST-II 
[20]. The predictions of our scenario are labeled in the same way as in Fig. 1. The two panels 
correspond to two different values of X, and we took Case 2 for the DM-quark couplings Cq. In 
the left panel, the band for / = 1 TeV (not shown) is very similar to the one for f = 500 GeV. 



exclusion limits of XENONIOO, that is presently unexplained. 

In Fig. 4 we compare a^i in our scenario with the DM direct detection constraints, taking 
the same values of the parameters that we used for the relic density in Fig. 1. The calculation 
of (Jsi is a straightforward application of the standard formalism [21] to our composite scenario, 
that is defined by the effective interactions of rj with quarks and gluons. The detailed equations 
are reported in appendix C, and can be summarized by the following estimate, which shows 
explicitly the relative size of the different contributions: 



asi ~ 3.5 ■ 10 



A 1 



-^Ocm^ 



10GeV\^ /125GeV^ ^ 



rrir, 



rrih 



2f_ 



^ + ^Re ^0.04c„ + 0.11q + 0.18c, + 0.22 ^ 



(5.1) 



The first term in square brackets comes from Higgs-exchange, while the second term comes 
from the direct coupling between the DM and the light quarks and gluons (through a loop of 
heavy quarks) in the nucleon. Contrary to the relic density, the DM direct detection bounds 
are sensitive also to the couplings between rj and the light quarks, mostly the strange. The 
Fig. 4 shows (Tsi for the Case 2 defined in section 2: we vary Q = = Cf, in the range 
defined by Eq. (2.9) and we take Cu,c,t = 0. Notice that the effect of the derivative coupling 
of Eq. (2.2) is very small at low momentum transfer and it has been neglected. 

For A = 10~^ (left panel of Fig. 4), the first term of Eq. (5.1) dominates and asi is of the 
same order of magnitude as in the non-composite case. By comparing with Fig. 1, one remarks 
that the correct relic density can be reproduced by m„ ~ 10 GeV, with a DM candidate that 
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lies close to the region preferred by DAMA, CoGeNT and CRESST-Il, or by ~ 80 GeV, 
in a region already excluded by XENONIOO, or by > 150 GeV, with no constraints from 
present direct detection bounds. For A = 10^'^ (right panel of Fig. 4), the difference w.r.t. 
the non-composite case is more important. The correct relic density is reproduced e.g. for 
/ = 1 TeV and ~ 100 GeV (see Fig. 1). This candidate is compatible with the present 
XENONIOO bound, and it is within the reach of near future measurements with improved 
sensitivity. For such small values of A, on the contrary, the non-composite case cannot be 
probed in DM direct detection experiments. 

We remark that the various contributions in Eq. (5.1) can have opposite sign and (par- 
tially) compensate each other, reducing asi. Above we assumed conservatively that all con- 
tributions are positive, but one cannot exclude a cancellation, e.g. switching to a negative 
value of A. Therefore the DM direct detection bound is quite sensitive to the details of 
the theoretical model, and for specific values of the parameters our scenario can avoid this 
constraint. 

Let us make now a few remarks about DM indirect detection experiments. Bounds from 
direct detection - as shown in Fig. 4 - allow for a light DM candidate with a mass < 10 
GeV. In this regime, the DM annihilation in the early Universe is entirely driven by the 
process rji] ^ bb, with a required cross section of the order of (crVrei) — 3 ■ 10^^^ cm'^s^^. 
This value weakly depends on the DM relative velocity, and therefore it remains constant 
from the freeze-out until today, resulting into a sizable antiproton flux on the earth, produced 
by the DM annihilation into quarks that hadronize. We estimated this flux using the tools 
in Ref. [22]. As already pointed out in Ref. [23], this antiproton flux can be larger than 
the one observed by the PAMELA experiment [24], resulting in a constraint on the DM 
mass that is complementary to the direct detection one. One should keep in mind, however, 
the large astrophysical uncertainty of these observations, mainly due to the propagation of 
charged particles in the Galaxy. A mass smaller than 10 GeV is generically disfavoured, 
with the exact bound depending on the adopted propagation model. In addition, the DM 
annihilations can also leave a trace in the cosmic microwave background power spectra, the 
size of the effect being proportional to {aVrei) /'fn-q- Assuming that (crVrei) is approximately 
velocity- independent, one finds that the Planck satellite should reach a sufficient sensitivity 
to observe such effect, as long as < 50 GeV [25]. 

We remark also that when the DM annihilation process occurs close to the Higgs resonance, 
the Breit-Wigner enhancement mechanism is operative [26, 27], resulting in a large boost of 
the annihilation cross section today, w.r.t. its freeze-out value. This boost factor can be 
invoked to explain some recent anomalies reported in cosmic ray data [26, 28]. Notice that 
the same resonant enhancement may play a crucial role for the pair annihilation of DM 
particles into two, monochromatic high-energy gamma rays. This channel is particularly 
relevant, given its peculiar spectral structure that cannot be confused with the continuum 
astrophysical backgrounds. A value of the DM mass too close to the Higgs resonance can 
exceed the upper limit ((crv^g;)^^ ~ 4 ■ 10~^^cm^s~^ for ~ 60 GeV) recently established by 
the Fermi-LAT collaboration [29]. 
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Finally, let us compare the DM direct detection bounds with those coming from collider 
searches of monojet plus missing transverse energy. Indeed, the singlet t] might be pair- 
produced through its coupling to quarks; while the two singlets escape the detector, these 
events can be identified by the initial state radiation of a gluon or a photon. In some cases 
[30] these searches can provide comparable bounds to DM direct detection experiments, in 
particular for light fermionic DM (with ^ 10 GeV). This is however not the case for 

our composite singlet for three reasons. First, the reach of monojet searches is reduced for a 
light scalar DM than it is for a fermionic DM [31]. Furthermore, in our scenario, the effective 
coupling of rj to gluons depends on the momentum of the process (at the momenta involved in 
DM direct detection experiments, t] couples more strongly to gluons than at collider energies), 
contrary to Ref. [30] where the effective coupling was assumed to be constant. Finally, in our 
case bounds on the invisible Higgs decay (section 4) provide the strongest collider constraint. 
For these reasons we neglect bounds from monojet searches in what follows. 

6 Combined results 

We present in this section a combined analysis of the parameter space for the composite sin- 
glet DM candidate rj. The results are displayed in Figs. 5-7, in the plane (m^. A). We require 
to reproduce the observed value of the DM relic density (the purple contour corresponds to 
fl^i = ^dm), and we take into account DM direct detection experiments (the green region is 
disfavoured by XENONIOO) and Higgs searches (the red region is disfavoured by ATLAS and 
CMS). 

• mh = 125 GeV, / = 500 GeV: a DM candidate with ~ 70 GeV 

We begin by taking the Higgs mass value presently preferred at the LHC, nih = 125 GeV, 
and choosing the characteristic pNGB scale close to the lower bound coming from electroweak 
precision tests, / = 500 GeV. We focus on Case 2 (a vanishing coupling q between the top 
quark and the DM), since it allows for a light r], < rrih. This scenario is illustrated in 
Fig. 5. 

In the non-composite case (/ = oo, dashed contours), a viable DM candidate lies in 
correspondence to the Higgs resonance (m,, — 60 GeV, A < 0.02), or beyond the kinematical 
threshold for annihilation into electroweak gauge bosons (m,, > 80 GeV, A ~ 0.04). 

In the composite case the situation changes considerably. Note first that the region 
^ ?7i/i/2 is disfavoured by the LHC Higgs signal excess, because the Higgs decay width is 
dominated by the hidden decay channel into DM, except for A ~ m\/{2p) ~ 0.03, where a 
cancellation in takes place. Note also that the direct detection bound is strengthened by 
the composite interactions, because of the contribution of direct DM— quark couphngs, see 
Eq. (5.1). 

Taking into account these constraints, a definite prediction follows from Fig. 5: ^ 70 
GeV with A < 0.02. At small values of A, the rj-rj-h vertex is dominated by the derivative 
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Figure 5: The contour = Qdai (solid dark purple line) in the plane {nij^, X), for = 125 
GeV, f = 500 GeV, assuming Case 2 with Cb = 1/2. The green shaded region is disfavoured 
by XENONIOO, the region delimited by a blue line is favoured by DAMA/CoGeNT/CRESST- 
II, and the red shaded region is disfavoured by the Higgs signal at the LHG. The solid light 
purple/ green/blue lines correspond to the same observables for maximal Cb (a = b = 1 in 
Eq. (2.9)). The dashed purple/ green /blue /red lines correspond to the same observables in the 
non-composite case, f = oo. Finally, the region below the yellow dot-dashed line corresponds 
to the theoretical preferred region defined by Eq. (2.6). 



coupling, therefore the rehc density becomes independent from A and the purple lines in 
Fig. 5 become vertical. For higher DM masses the derivative interaction becomes too strong 
to accommodate the relic density (unless one enters in the "cancellation region" , at relatively 
large values of A, which is however excluded by XENONIOO). 

One may ask if the region of parameters that is allowed phenomenologically is also com- 
patible with the theoretical expectations. Independently of the specific model, we expect from 
Eq. (2.6) that > A/^. The region satisfying this relation lies below the yellow dot-dashed 
line in Fig. 5, and it is compatible with the phenomenologically preferred region. 

• mh = 125 GeV, / = 1 TeV: DM candidates with ~ 60 GeV and 100 < < 500 GeV 
As the scale / increases, the composite interactions become weaker, and the bounds from 
the LHG Higgs signal and from XENONIOO become less stringent and closer to the non- 
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mi, = 125 CeV, f = l TeV. case 1 



m,, = 125 GeV, f = 1 TeV, case 2 




nir, [GeV] 



[GeV] 



Figure 6: The same as in Fig. 5, but with / = 1 TeV and with a comparison between two 
scenarios for the top quark couplings: Case 1 (left panel) and Case 2 (right panel), as defined 
at the end of section 2. 



composite case. This is illustrated in Fig. 6, where we take / = 1 TeV. In particular, all 
values of are viable for A < 10~^. 

The correct DM relic density can be accommodated for lying a bit below or above the 
Higgs resonance at m/i/2 ~ 60 GeV, or for > 100 GeV, where the derivative interaction 
Tj-rj-h in Eq. (2.2) becomes of the right order to give the correct annihilation cross-section 
above the WW threshold. Furthermore, for relatively large values of A, one enters in the 
cancellation region described in section 3: the DM annihilation is suppressed and the relic 
density can be accommodated even for very large values of the DM mass, up to ^ 500 
GeV in Case 2 (right panel of Fig. 6). If the annihilation into ti is stronger (Case 1, left panel 
of Fig. 6), the allowed region closes earlier, at rrir, — 200 GeV. 

As discussed above, composite models prefer A < rn"^/ f^ (the region below the yellow 
dot-dashed line) that is compatible with the Higgs-resonance region for A < 0.003, and with 
the region dominated by the derivative coupling, for A < 0.02. On the contrary, the cancel- 
lation region is slightly disfavoured theoretically, even though A larger by a factor of a few is 
sufficient to realize the cancellation. 

• mh = 145 GeV, / = 500 GeV: DM candidate with < 10 - 20 GeV and ~ 60 GeV 

In case the LHC excess at 125 GeV were not confirmed, the Higgs boson might be heavier, 
as long as it decays invisibly with a sufficient rate to avoid the LHC bound. In order to 
illustrate this possibility, in Fig. 7 we choose a representative value nih = 145 GeV, assuming 
for definiteness / = 500 GeV and Case 2 (the results are very similar in Case 1). 
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Figure 7: The same as in Fig. 5, hut with nih = 145 GeV. 



The LHC bound is satisfied easily below the kinematical threshold for Higgs decays into 
DM, < m/j/2, in a region of parameters which is very much complementary to the one 
allowed in Fig. 5. As a consequence, a light DM candidate with a mass < 10 GeV is 
compatible with XENONIOO and LHC. A light singlet with 10 GeV < < 20 GeV could 
in principle explain the DAMA/CoGeNT/CRESST-lI results (of course, in tension with the 
XENONIOO bound). However the value A < 0.3 required by the relic density is slightly 
smaller than the one needed to fit the signal in these experiments: for larger A, r] accounts 
only for part of the DM relic density. 

In any case, the large coupling A > 0.1, needed to explain the relic density when < 20 
GeV, is in contradiction with theoretical expectations. Indeed, the singlet receives at least 
a contribution to its mass-squared of order Af^ > (80GeV)^, and the NDA expectation is 
actually much larger, > \p > (160 GeV)^. Therefore, the case of a light singlet is 
disfavoured in the context of composite models, since it requires a large cancellation between 
different contributions to the singlet mass. 

The other viable DM candidate in Fig. 7 lies just before the Higgs resonance, at ~ 60 
GeV, as long as A < 0.02. This solution lies in the theoretically favoured region, and it 
requires a not too large R.e{cb), in order to avoid the XENONIOO bound. 

In summary, it is interesting that the composite scalar DM model is very predictive, despite 
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the new non-renormalizable interactions, generated by the strong sector at the TeV scale. In 
some cases, the composite singlet DM is even more constrained than the non-composite one. 
Indeed, in good approximation, the relic density and the LHC bounds depend only on m,,, A 
and the derivative interaction between rj and the Higgs, whose strength is uniquely fixed by 
the pNGB scale /. Therefore, our results are independent from the specific composite model, 
up to the relatively small modifications due to the choice of the parameters q (compare the 
two panels in Fig. 6) and Cb (compare the dark and light contours in Figs. 5-7). The largest 
model-dependence lies in the DM direct detection bound for small values of A. The most 
interesting DM candidates (for rrih = 125 GeV) have mass ~ 60 — 70 GeV, around the 
Higgs resonance, or lie in the region 100 GeV < < 200 GeV. In both cases small values of 
the Higgs portal coupling, A < 0.03, are favored both phenomenologically and theoretically. 

7 Conclusions and prospects 

We have studied the phenomenology of a composite DM candidate, arising from the 0(6) /0(5) 
global symmetry breaking pattern. The light scalar spectrum of this model is formed by five 
pNGBs, which correspond to the Higgs doublet H and the real singlet t], that plays the role of 
the DM. They are lighter than the dynamically-generated scale rup ~ TeV, due to their pNGB 
nature. We have shown that the DM phenomenology is strongly affected by the composite 
nature of the pNGBs, that allows for 

(i) derivative interactions between rj and H which grow as / P ■, where p is the relevant 
momentum; 

(ii) contact interactions between r] and the SM fermions ip that break the ?7-shift symmetry 
and are suppressed by m^p/p; 

(iii) modification of the Higgs couplings to gauge bosons and fermions at order v'^/p. 

We found that this scenario is very predictive, and it is actually easier to test than the 
non-composite case, which depends only on the Higgs portal interaction A. 

The most important phenomenological consequences of the properties (i)-(iii) are the 
following. The i] annihilation cross section grows with because of the derivative interaction 
between H and 1]. Therefore, to reproduce the DM relic density, one needs m,, < 500 GeV 
(for / < 1 TeV). The same derivative interaction implies that the invisible Higgs decay in two 
DM particles is possible even for vanishing A. On the other hand, the composite Higgs can 
be slightly fermiophobic and privilege decays into gauge bosons, as preferred by recent LHC 
data. The couplings can significantly increase the DM scattering cross-section on nuclei: a 
signal in DM direct detection experiments is possible even for vanishing A, as long as 10 GeV 
<m„< 100 GeV and / < 1 TeV. 

Once all the constraints are taken into account, we are left with only three regions of 
parameters for our DM candidate: 
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1. Light mass region, with ~ 10 — 20 GeV, close to the values preferred by the signals of 
DAMA/CoGeNT/CRESST-II. In this case the Higgs has a substantial invisible width 
into DM, that would be strongly disfavoured if the LHC hint for rrih — 125 GeV were 
confirmed. Such light DM also contrasts with the theoretical expectations, because it 
requires a large coupling A > 0.1, and therefore a strong cancellation is needed to make 
the DM mass small. 

2. Intermediate mass region, with ~ 50 — 70 GeV, where t] annihilates mostly through 
the Higgs resonance. This region of parameters is almost the same as in the non- 
composite case, but a discrimination is possible for / ~ 500 GeV, where the LHC 
and DM direct detection constraints are significantly different. For < m/i/2, Higgs 
invisible decays are expected, resulting in a reduction of the LHC Higgs signal. 

3. Heavy mass region, with ~ 100 — 500 GeV, in which the annihilation cross-section 
is dominated by the derivative interaction between H and t]. The latter is fixed by the 
value of /, and the relic density is actually too small for / ~ 500 GeV, but it becomes 
the correct one when / ~ 1 TeV. This case is possibly the most interesting, since the 
relevant DM interactions are completely determined by its composite nature, and it can 
be tested in near future DM direct searches (in contrast to the non-composite case). 
Moreover a DM mass close to the Higgs mass would strongly suggest a common origin, 
that is to say, both are generated radiatively by order one couplings to the top quark. 

Other indirect signals of our scenario are the following. Since the Higgs is also composite, 
its couplings to the SM particle differ from an ordinary Higgs, as it can be explicitly seen 
in Eqs. (4.3)-(4.7). It has been shown in Refs. [8, 32] that the LHC with greater integrated 
luminosity could discriminate between standard and composite Higgs for values of / ~ 500 
GeV. Furthermore, the Higgs can decay into two r] that could be measured as an invisible 
decay width, as long as t] is sufficiently lighter than the Higgs. When t] is heavier, at large 
center-of-mass energy searches for monojets plus missing transverse energy might be the 
most sensitive to the derivative coupling, and they deserve further investigation. Also, heavy 
resonances p of the strong sector could decay into the DM with a sizable branching ratio for 
large Qp. For example, the decay p — )■ rjriW could result in a distinctive signal with leptons 
and large missing energy in the final state. 

Let us briefly comment on another cosmological implication of a singlet scalar. It has 
been pointed out [33] that a real singlet with the parity of Eq. (2.1) can help in inducing 
a strongly first order electroweak phase transition, and thus play an important role in elec- 
troweak baryogenesis, as well as leave an observable spectrum of gravitational waves [34]. It 
would be interesting to check whether this possibility can be realized together with the singlet 
being DM. Although it seems unrealistic for an elementary singlet (the couplings needed to fit 
the relic density are too small to play a role during the phase transition) , it could be possible 
for a composite 7], that has large additional interactions. 

In summary, we have presented a minimal composite framework that accounts for the 
Higgs boson and for a scalar DM in a natural way. Future DM searches (e.g. XENONIT) and 
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LHC measurements shall be able to close on the three most promising regions of parameters 
described above, and thus establish or refute the connection between the EWSB sector and 
the DM sector. 
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A The composite 0(6)/0(5) model for dark matter 

In this appendix we present a model of a composite sector that accommodates the Higgs 
doublet H and the DM singlet r] as pNGBs. The analysis is intended to provide the reader with 
a concrete realization of the effective lagrangian that we used all over our phenomenological 
study, and to compute the expected values of the pNGBs masses and couphngs. The coset 
S0{&) / SO{b) provides the minimal realization of a set of pNGBs with the quantum numbers 
of H and rj. This coset was first studied in Ref. [3], with focus on the modifications of the 
Higgs phenomenology in the presence of a real scalar singlet. Here we demand that t] plays 
the role of DM particle and we derive its properties, under the requirement of a consistent 
EWSB. 

We assume a strong sector, with a mass gap mp ~ TeV, whose global symmetry breaking 
pattern is SO{Q) 50(5). This can be achieved by a composite field S in the fundamental 
representation 6 of 50(6), that acquires a VEV 

So = (0,0, 0,0, 0,1)^. (A.l) 

The five NGBs transform as a 5 of 50(5), which decomposes as 4 © 1 ^ (2, 2) © (1, 1) under 
the subgroup 50(4) ~ SU{2)l x SU{2)r. The NGBs describe the fluctuations along the 
broken directions, whose generators can be written as 

Ti = -^{SfS^-Sp!), a = l,...,5. (A.2) 

The 50(6)/50(5) coset is parametrized by 
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with tt" = {hi, 7]), n = yJ2i=i ^1 + ^^"^ ^ ^^1^ redefinition 

, hi n f] TT 

hi = —sm- , ?7 = - sm - , (A.4) 

with /i^ = /i^. The usual Higgs doublet is given by if = f[{hi + ih2)/V^, {h3 + ih4)/\/2]. 
In the unitary gauge three NGBs are eaten by the weak gauge bosons, and one is left with 
the Higgs boson h = h^ and t]. 

The NGB chiral lagrangian reads, up to terms with four derivatives, 

^,2 _ .^2 , ^^2 , {hd,h + r^d,vr 



Q^h^ ( 1 



^ 2cos2^ 



(A.5) 



From the last term in Eq. (A.5), it is clear that the scale of EWSB is set by the VEV of h: 
(h) = v/ f = a/^. Note that Ckin is symmetric under the parity : — > — r], but higher 
derivative terms in the SO{Q) / SO{b) chiral lagrangian, involving the Levi-Civita tensor (such 
as the Wess-Zumino-Witten term [3]), are not. Since is crucial to make r] stable and thus a 
viable DM candidate, we need to assume that it is a symmetry of the whole composite sector 
and it is not spontaneously broken. This amounts to take the global symmetry breaking 
pattern to be 0(6) 0(5), with 

P, = diag(l, 1,1, 1,-1,1), (A.6) 

and to require {rj) = 0. Note also that, once h acquires a VEV, its kinetic term receives a 
correction from the third term in the square bracket of Eq. (A.5). The physical Higgs and 
DM bosons (with a canonical kinetic term) are defined as follows: 



T'phys h (/l) ^phys 



V ■ (A. 7) 



In all the observables we studied, we consistently used the couplings of the physical fields 
(but dropping the subscript 'phys^ everywhere). 

The EWSB is triggered by the couplings of the SM gauge bosons and fermions to the 
composite sector, that break 0(6) explicitly, generating an effective potential for the pNGBs 
at the one-loop level. Still, we will show that under some conditions Pr, is preserved by the 
gauge and fermion interactions as well. In this case, the most general form of the effective 
potential for h and 1] can be written as 

Such potential must be minimized taking into account the constraint K^+yf = sin((7r)//) < 1, 
that follows from Eq. (A.4). We look for a minimum with 0<h = v/f<l,to realize EWSB, 
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and with 77 = 0, to preserve DM stability. The necessary and sufficient conditions for this to be 
a local minimum are < —fil < pXh and fx'^Xh > /^^A, where we neglected the dimension-six 
terms in Eq. (A. 8).^ Then, the pNGB physical masses are given by 

ml ~ -2fil = 2XhV^ , ml ^ jj^ + Xv^ , (A.9) 

up to O(^) corrections. 

The size of the effective potential coefficients depends on the way the SM gauge bosons and 
fermions couple to the composite sector. Note that, in order to reproduce the hypercharge 
of the SM fermions, the global symmetry of the composite sector should be enlarged to 
0(6) X U{l)x, with the hypercharge defined by 

r = r| + x. (A.IO) 

Since 77 is a gauge singlet, its NGB nature is not affected by the gauging of SU{2)l x [/(l)y, 
that is to say, it does not acquire a potential through the gauge loops. Therefore, the elec- 
troweak gauge bosons only contribute to the Higgs potential, inducing /i^ ~ 9^f^pl (167r^) and 
Xh ~ 5'^/(167r^) (see Ref. [3] for details). These contributions are generally smaller than those 
coming from the top quark, which will drive EWSB. 

The Eqs. (A. 5), (A. 7) and (A. 8) define all the relevant interactions among 77, h and the 
gauge bosons. The corresponding Feynman rules are displayed in Table 2. 

Coming to the SM fermion interactions with the strong sector, we assume the partial 
compositeness scenario [35, 2], which is preferred by the constraints on flavour violation: each 
SM chiral fermion t/i couples linearly to a composite operator of the strong sector: 

Cint = X^ip + h.c. . (A. 11) 

This leads to a mixing of the SM fermions with the heavy composite states (of mass ~ mp) of 
order X^/gp. By repeating this reasoning for each SM chiral fermion, the Yukawa couplings 
turn out to be given by ~ X^^X^^^/gp. 

Following standard techniques [2, 37], we will promote the SM fields ip into spurious \E' 
forming complete SO (6) representations. The \E' are defined to transform as the 5*0 (6)- 
multiplet in such a way that the interaction Eq. (A. 11) can be written as an invariant 
under 5*0(6). The interactions between the SM fermions and the pNGBs h and r] will crucially 
depend on the representation of under 50(6). We will identify the simplest representations 
that (i) allow to generate the Yukawa coupling, (ii) induce a pNGB effective potential that 
realizes EWSB satisfactorily, (iii) preserve the parity that guarantees the stability of our 
DM candidate rj. 

In addition to these necessary requirements, the choice of the 5*0(6) representation for \E' 
also determines whether breaks or preserves the 77 shift symmetry U{l)rj, that is to say, 
whether or not V^-loops contribute to the effective potential for rj. This is crucial to determine 

The requirement for this to be the global minimum (in the region + rj^ < 1) involves some extra 
lengthy conditions, that are satisfied in a wide range of the potential parameters. 
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Table 1: 5*0(6) representations for the spurions \E' that can embed the SM fermions ip preserv- 
ing P^. The symbol y/ (x) indicates that the Pj^-preserving embedding does not break (does 
break) the 5*0(2)^ subgroup. 



the DM mass. To study this issue, it is useful to decompose the SO{Q) multiplets under the 
maximal subgroup 50(4) x ^0(2)^^ ~ SU{2)l x SU{2)r x f/(l)^, where ^0(2)^ ~ ?7(1)^ is 
precisely the symmetry associated with the NGB r/, that is generated by = defined in 
Eq. (A.2): 

T, = i= f \ . (A.12) 



V2 V 02x4 0"2 

From Eq. (A. 6) and Eq. (A.12) we have 

[P„T,]^0, (A.13) 

therefore one cannot assign to the SM fields a definite parity and a non-zero U{l)r, charge 
at the same time. This means that both symmetries can be preserved by the SM couplings 
to the strong sector only if the SM fields transform trivially under the 5*0(2)^, i.e. they are 
not charged under f/(l)r,. 

Let us classify the spurions ^ that can accommodate the SM fermions while preserving P^. 
The SM fermion isosinglets (isodoublets) can be embedded in any 5*0(6) representation that 
contains a singlet (doublet) of SU{2)l. We find that an embedding preserving P^ is possible 
using a vector representation, \1' ~ 1, 6, 15, 20', . . . , or a pair of spinor representations, 
4 + 4, 10 + 10, and so on (when acting on spinors, P^ interchanges conjugate representations). 
Their SU{2)l x SU{2)r x f/(l)^ decomposition reads 



4 = (2,l)+i©(l,2)_i 

6 = (2,2)o©(l,l)+2©(l,l)-2 , 

10 = (2,2)o©(3,l)+2©(l,3)-2 , 

15 = (l,3)o©(3,l)o©(l,l)o©(2,2)+2©(2,2)_2 , 

20' = (3,3)o©(2,2)+2©(2,2)_2©(l,l)+4©(l,l)-4©(l,l)o • (A.14) 

The components that can contain the SM fermions and do not transform under U{l)n are, 
for the isodoublets, the (2, 2)o in the 6 or in the 10, and for the isosinglets either the (1, l)o 
or the (1, 3)o, which are present in the 1, the 15 or the 20'. These results are summarized in 
Table 1. 

In the following we will describe the properties of two models. In the first one we will embed 
all the SM fields in the 6 of 50(6), and therefore the right-handed couplings to the strong 
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sector will break the U{l)r^ symmetry explicitly. The top quark loops will give the largest 
contributions to all terms of the pNGB effective potential in Eq. (A. 8); this corresponds to 
Case 1 of section 2. In the second model, on the other hand, we will embed the up-type 
quark singlets, and in particular tn, in the 15, so that U{l)r, will be preserved by the top 
couplings; doublets and down-type quark singlets will be embedded in the 6 as before. In 
this case the DM particle can be lighter than the EW scale leading to a substantially different 
phenomenology; this corresponds to Case 2. Both models satisfy the properties (i)-(iii) given 
above. 

Before going to these two specific models, let us briefly mention an alternative way to keep 
T] light. If the SM embedding in S0{6) representations preserves the S0{5) subgroup, that is 
broken by the Higgs VEV down to an 50(4)', we would have after EWSB four exact NGBs, 
the three SM ones eaten by the weak gauge bosons plus t]. In this case the mass of rj would be 
protected by this 50(4)' symmetry and not by the 5*0(2)^ of Eq. (A. 12). It is clear, however, 
that the SM fermion doublets cannot be embedded in complete 50(5) representations, so t] 
gets always a mass from loops of this sector. A relative light t] could be achieved by requiring 
the smallest possible coupling of the fermion doublets to the strong sector. Here we do not 
explore this possibility any further. 

A.l Case 1: all the SM fermions in the representation 6 

Let us begin with embedding all SM fermions into \Ef ~ 6 [3]. A SM electroweak doublet can 
be embedded in the the bi-doublet with = +1, while for a SM isosinglet the only embedding 
with = +1 (to avoid the breaking of P^) is given by = (0, 0, 0, 0, 0, ipR)- Coming to the 
U{l)x assignments, recall that the pNGB field E does not carry X-charge. Using Eq. (A. 10), 
the isosinglet up-type and down-type quarks should respectively be embedded into a spurion 
with = 2/3 and = —1/3. Then, in order to allow for the Yukawa couplings to 
E, the isodoublet quarks should be embedded into two spurious, one with = 2/3 to 
generate the up-quark Yukawa, y„ ~ \qj^\up,/gp, and one with Xqi^ = —1/3 to generate the 
down-quark Yukawa, i/d — Xqi^Xd^^/ Qp-^ Here and in the following we tacitly assume that the 
lepton couphngs to the composite sector are generated in analogy with the quark ones. 

As can be seen from the U{l)r^ charges of the 6 components in Eq. (A. 14), the rj shift 
symmetry is broken by the embedding of the right-handed quarks. Therefore rj will receive 
an effective potential from loops involving only isosinglet quarks. This can be seen explicitly 
by writing all the invariants that involve the quarks and the pNGBs. The techniques of 
Refs. [36, 37] are useful to count the number of possible invariants, that we write using the 

^ In the case of the third family, the coupUng A^^ imphes a correction to the SM value of the coupling 
ZbLbL, because it corresponds to T^{bL) = 1/2 different from TKb^) = —1/2 [38]; we assume a small value 
Xq'^ ^ yb, to guarantee a small enough correction to Zb^bL. 
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j)'W(g'B 'interactions 
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n ^ 
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-2ivXVl-^ iv(k,+ k,)2 


-2iv A, 


\ /■'Pi 
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-2il(l- U(Pi+ Pa)' _ p,. p, 
f (1- s ) f (1- ? ) 


-2il 




-3im^ c 2iv TCd + 0^)^-0 d 1 


V 


J}14i^'B/S'M Int. 


Composite 


ifon-Composite 


n ^ / 77) 


V^ivy^ [Re (c^) + i Im (c^) ] 0^ e 6 





:< 


t/^j,e15 








•< 






■< 




i^triw V = W 




^.^^zVl-^ V = Z 


^d^z V = Z 



Table 2: Feynman rules for the relevant interactions between the Higgs boson h and the 
DM particle t] (upper panel) and for the interactions of h and rj with the SM fermions and 
gauge bosons (lower panel) . The second column corresponds to the composite model with coset 
0{6)/0{5) , while the third column corresponds to the SM plus a real scalar singlet. 
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spurion notation 



J = UR ^URil -rf- , 

(%S)^(E^v|/^, ) = ]^dL]^dLh\ (A.15) 

We can now justify the parameterization of Eq. (2.9) for the DM— bottom couphng. The first 
and second hnes in Eq. (A.15) generate the quark masses and, at order rf, they also give a 
contribution to the couphngs of Eq. (2.3), Ac^ = Aq = 1/2, up to C(0 corrections. The 
third and fourth hnes in Eq. (A.15) correspond to derivative couphngs analog to the one in 
Eq. (2.5), that can be cast in the form of Eq. (2.3) by using the Dirac equation. Using NDA to 
estimate the coefficients of these invariants, one obtains /S.Cu,d = 0{\l^^^^^/gp)+i 0{\^^^^^/ gj,). 
Note that the ratio ^s^^^dnldl varies approximately in the range \y\^Jg'],-, 1], depending on the 
relative size of \i^,q'^ and Xur^r- Adding the two contributions to the DM— fermion couplings, 
one recovers Eq. (2.9).*^ 

The Feynman rules for the interactions of t] and h with fermions, following from Eq. (A.15), 
are reported in Table 2. 

The effective lagrangian in momentum space, obtained after integrating out the composite 
sector, can be written as 



- [M"(p2)(^,,S)(S^vl/„J + M<^(p2)(^r^, S)(S^xl/,J + h.c] , (A.16) 

where p is the momentum of the fermion fields. Parametrically, at small momentum {p < rrip) 
we have 

n^-i, n^~4' M"~^^^/, M'^-^^/. (A.17) 

9p 9p 9p 

Eqs. (A.15) and (A.16) imply that the fermion masses read rriu^d = M"''^(0) / \/2(v/ f)\^l — ^. 
Loops of SM fermions contribute to the Coleman- Weinberg potential for the pNGBs. Since 
the heaviest fermions give the dominant contribution, we concentrate on the third family of 



^ In general, the contributions of order A^^ dn/dp could have negative sign and partially compensate the 
contribution 1/2. In our analysis we did not consider the possibility of such a cancellation. 
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quarks. At one-loop, from Eq. (A. 16) one obtains 

v,f{h, rt) = -2Nc / (Ifi {log [p' + nf - ^h^) (n*« - nl-(i -r,^- h-)) 

+ log[(gi,tfi)o(gl,6H)]| , (A.18) 



|MT,o, O ,0 



2 

where the overall negative sign comes from the fermion loop, the factor of 2 from the fermion 
polarizations and Nc is the number of colors. 

For a momentum larger than the mass of the resonances mp, one expects the form factors 
and M* to fall off fast enough to make the integral in Eq. (A.18) convergent. At p < rup, 
for <^ Qp one can expand the logarithm in Eq. (A.18), since each power of /i^ and rf is 
associated with a suppression factor (A^A^'/^'p)- Oiis can easily check that all terms in the 
expansions in h"^ and rf receive contributions from the top couplings Ag^ and At^, therefore 
one can safely neglect the terms in A^^ and Af,^, that we assume to be of the order of the 
bottom Yukawa. The first logarithm in Eq. (A.18) provides the dominant contribution, coming 
from the top-quark, to the effective potential K// defined in Eq. (A. 8). By expanding this 
contribution and integrating over momenta, we obtain 



2 Ncml . ^ ,2 N 



IGtt 

Nr /_ (2)^4 

R 



\ - 1^(2^2 A, 



A 



^,i2c,y^9l + 2c?hl) , (A.19) 

where q and c\ are order one coefficients that account for the uncertainty in the integrals over 
the form factors. More specifically, the q come from the second order in the log expansion, 
and therefore they involve an integral over a higher power of the form factors, that in some 
explicit models is suppressed (see e.g. Ref. [39]); in the following we neglect for simplicity 
these terms, when compared to the q ones. 

In order for C, = —fi\/{XhP) to be smaller than one, a partial cancellation is needed 
between the terms in that requires A^^ ~ A^^.' Then, from Eq. (A.19), using yt ~ 
Xq^Xtj^/gp and Eq. (A. 9), the Higgs mass can be written as 

ml = e ml ^ C3(250 GeV)^ (^)' . (A.20) 

The coupling A is generated by the same term that controls also Xh, so we expect it to be 
positive and of the same order 



^ We are neglecting gauge contributions that are smaller than the top one. 
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(for nih = 125 GeV this corresponds to A ~ 1/16). The singlet mass can be written as 

,2 ,.2 „ .2 _2^_2C2'^? 



Since the last factor in Eq. (A. 22) is larger than one, rj is generically heavier than the Higgs. 

A. 2 Case 2: the right-handed top quark in the representation 15 

As mentioned above, if t/j is embedded into an SO{Q) representation preserving ?7(l)r,, the 
DM 7] is naturally lighter than the Higgs. In this section we provide an example of such a 
model, where up-type isosinglet quarks are embedded in the 15.^ As in the previous section 
we keep all isodoublets and down-type isosinglets in the 6. The up-type Yukawa can be 
generated when ur is embedded in the (l,3)o C 15. Taking for definiteness Xq^ = 2/3, in 
the up-type quark sector Eq. (A. 15) is replaced by 



(fr^,S)#(S^M/,J = -UL^ULh'. (A.23) 

As expected, no terms involving t] are generated, therefore the DM— top coupling q of Eq. (2.3) 
vanishes, and there are no loop contributions to the effective potential for rj from the up-quark 
sector. 

There are still, however, DM couplings to the down-type quarks; the heaviest is the 
bottom, therefore Eq. (A. 22) is replaced by 

"^?^/^? = ^2^^'X' (A.24) 

which results in a DM mass > 0(10) GeV (Aj,^ > ?/f, — 1/40). Notice that, by embedding 
also the Br in the 15, we could make the DM even lighter. Note also that the DM— Higgs 
coupling A is correlated to the DM mass, 

in agreement with the general expectation of Eq. (2.6). Therefore, the lightest the DM is, the 
weakest its Higgs portal interaction becomes. 



^ Embedding ta in the 1 of SO{6) is the simplest possibihty for hght DM, but it is not satisfactory for 
EWSB, because in this case Ir does not break the Higgs shift symmetry either. As a consequence, one needs 
to invoke a cancellation in the strong sector to obtain w ^ /, and moreover the Higgs mass tends to be too 
small. 
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Similarly to the previous model, starting from Eq. (A. 23) one can compute the top quark 
contribution to the effective potential, 



(2^ 



log 



P 



■'■■'^0 



TT'' 



TT* 



\M 



t\2 



(A.26) 



Expanding and integrating the logarithm, one finds 

i2 



^ p 
167r2 

167r2 



'^1 ■^QL + 0*^2 \ii + Cg lO, 



(2)x4 



.(2).,4- 



mr 



(A.27) 



where q are coefficients of order one. As in the previous model, in order to obtain f ^ / 
we need a partial cancellation between the different terms of fil, that can be achieved if 
\l ~ -^tfl ~ \fyWp (and if Ci and C2 have opposite sign). In this case the dominant terms in 

(2) 

the quartic Higgs couplings are expected to be those proportional to c]^ 2 ~ 1) that lead to a 
Higgs mass 

{A.28) 



2 -^c 2/- 2 ^1'7nr1■^A^2 ^ / 



4,4 1/.' e = (170 GeV)^ 



0.1 V2 TeV/ ■ 

In summary, EWSB can be smoothly obtained when is embedded in the 15, with a Higgs 
naturally close to the presently favoured low-mass region. 



B Dark matter annihilation cross sections 

In this appendix we outline the calculation of the relic density and collect for completeness 
all relevant cross sections. Following the standard recipe [9], the central ingredient in the 
calculation of the DM relic density is the thermal-average of its total annihilation cross section 
times the relative velocity Vrei, 

m f°° I 

{avrei) = ^ / rfsa(s)Ws - 4m2/Ci (xA/s/m^) , (B.l) 

047r Xn^Q Jirn^ V 

where a{s) = sVreiJ^F^ivn ~^ -^)) being a{riri — )• F) the annihilation cross section into the 
SM final state F. In Eq. (B.l) the equilibrium number density for r] is given by 

"^^^2(x) 

Here /Ci,2 are the modified Bessel function and x = rrirj/T, where T is the temperature. A 
Boltzmann equation describes the evolution of the number density n{x) of the DM particle 
during the expansion of the Universe, 

f = -^(--)(Y^-Yy , (B.3) 
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where the yield is defined by Y = n(x)/s(x), and the entropy density s{x) is given by 



Six] 



45x3 



(B.4) 



The Hubble parameter evaluated at x = 1 is given by 



H 



45 MpL 



(B.5) 



the value of the Planck mass is MpL = 1.22-10^^ GeV, and is the effective number of degrees 
of freedom. The integration of Eq. (B.3), from x = to the present value Xq = m^/Tg, gives 
the DM yield today, Yq, which in turns is related to the DM relic density [10], 



, _ 2.74 ■ lO^m^Yp 
" ~ GeV 



0.1126 ±0.0036 . 



(B.6) 



The computation of the annihilation cross sections makes use of the Feynman rules col- 
lected in Table 2. The relevant final states are the heavy fermions (the t, 6, c quarks, as well 
as the r lepton), the electroweak gauge bosons and the Higgs. Defining = avreio^iw ~^ 

F), we find 



y/l-4m2/s 



X 



Re(c^)^ + 



slm(c^)^ + {s - Ami 



r(l-20^ 



V'JJ615 

ar,r,{hh) 



4(1 -0^ 
N.mlis - 4m2 )3/2r 
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(B.7) 
(B.8) 

(B.9) 

(B.IO) 



1 — Am\/s 



X 



16/*7r(l — i'^Y{s — m\) \m\ + m^(s — Am\ 

^Pi + - 4^«h)] arctanh (A) [ml - 2Xp{l - 

(s - 2m2)2A 



A 



where we defined the ratios 



2mi 



{s - mlY + Tlml 



(b'.ii) 

(B.12) 
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The lengthy expressions for the functions and B^^ are hsted at the end of this appendix. 

Notice that taking the hmit / — )■ oo, we recover the known expressions for the annihilation 
cross sections in the non-composite case, 
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(B.13) 
,(B.14) 
(B.15) 
(B.16) 

(B.17) 
(B.18) 
(B.19) 



function of the DM mass, for 



We plot the annihilation cross sections in Fig. 
Vrei = 2/3, A = 10~^ and ruh = 125 GeV, taking / = 500 GeV in the left-hand panel and 
/ = cxD in the right-hand panel. For comparison we also show (black dashed line) the typical 
value for (aVrei) suggested by Eq. (3.1), needed for the thermal freeze-out of the DM particle. 
In the left panel of Fig. 8 one can distinguish the four DM mass regions described in section 
3, that characterize the qualitative behavior of the composite singlet relic density. 

Finally we report the expressions for the coefficients and Br^ri in Eq- (B.ll): 



Ar,r, = f{ml[AfX{^-im-2) + i5^-A)s]-mU[s + 2fXil~0\ 
+ (1 - Os [^f'X'il - - 2fX{l - 20s + s'] + 3m^e } , 



B. 



vv 



mr 



n=0 



(B.20) 
(B.21) 
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m, [GeV] m, [GeV] 

Figure 8: Annihilation cross sections times the relative velocity Vrei for the DM candidate rj. 
We plot Eqs. (B.7)-(B.ll) for f = 500 GeV (left-hand panel), and Eqs. (B.13)-(B.16) for the 
non- composite case (right-hand panel); in both cases A = 10~^, rrih = 125 GeV, Vj-ei = 2/3. 
We vary the coefficient Cb as in Fig. 1. For the top quark channel, we show both Gase 1 and 
2 (dotted and solid cyan lines, respectively). The black dashed line indicates the benchmark 
value for DM freeze-out, (avrei) = 3 ■ 10~^^ cm^s"^. 



with 

90 = (l-OV{32/«A^(l-OT + m;.[2f A(2e-1) + .]'} , (B.22) 
= 2si^-l)^{my[2fXi2^-l) + s]+32f<^\%^-in[s-fXi^-l)]}, 

(B.23) 

54 = (e - 1)' {32f A^(e - 1)^ - 128 f\%^ - l)es + ifX'lim - 1) + 1]^^'+ 
AfX{2^ - l)s' + - m^s {A8fX' (2^^ _ 3^ + i)^ + 

48/2A(e - 1)2(2^ - l)s + [^(11^ - 24) + 12]s^ } , (B.24) 
g, = 2i^-l){32fX'i^-lfe-8f^X\^-l)[2a^ + 2)-l]s~ 
4m2 [2fXi2^ -l) + s] [AfXi^ - - 1) + (5^ " 2)s] + 
2/2A[e(14e - 13) +A]s' + (3e - 2)s'} , (B.25) 

gs = 16/2A(e-l){/^A(e-l)[e(7e-4) + l]+4mJ(l-20e} + 

As {AfX{^ - l)[2(e - 2)e + 1] + mliS - 1^} + (3^ - 2)(5^ - 2).^ , (B.26) 

510 = 8e[2/2A(e-l)(3e-l)-2m;e+(e-l)s] , (B.27) 
912 = ■ (B.28) 
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C Dark matter scattering cross section on nuclei 



The spin-independent elastic scattering cross section of a real scalar t] on a nucleus can be 
parameterized as 

If m^m„ y [Zfp + {A- Z)fn]^ 

where is the neutron mass, Z and A — Z are the number of protons and neutrons in the 
nucleus, and the function fp (/„) describes the coupling between i] and protons (neutrons): 

« ^ m„ 27 m„ 

q=u,d,s ^ q=c,b,t ^ 

The hadron matrix elements /|^'^'* parametrize the quark content of the nucleon. They are 
extracted from pion-nucleon diffusion measurements. Following a recent analysis [40] we take 

/(f) = 0.017, 4^^ = 0.022, /(f = 0.053, 
= 0.011 , = 0.034 , fff = 0.053 . 

The DM couples also to the gluons in the nucleus, through a loop of heavy quarks (those with 
TTiq ^ Aqcd)- The corresponding matrix element is given by 

q=u,d,s 

In Eq. (C.2), the parameter describes the short-distance effective interactions between r] 
and a given quark, normalized as 

aqV^qq ■ (C.5) 

q 

In our composite model there are two sources for this effective interaction: the direct 
coupling Cg of rj to quarks, and a coupling mediated by the exchange of the Higgs boson. More 
precisely, only Re(cq) enters in a^, and thus contributes to cxsi, while Im(cg) leads to a coupling 
rfq^^q, that contributes only to the spin-dependent cross section. Coming to the Higgs in the 
t-channel, one should notice that the momentum transfer in the elastic scattering is very 
small, therefore the Higgs-exchange reduces to a contact interaction suppressed by m\. In 
addition, the derivative coupling between rj and h can be neglected, because it is proportional 
to g^, and one is left with the coupling A only. All in all, we find 
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Sl9i?6l5 - — ■Z:j2 • (C-"^) 

Putting all the ingredients togetb 
Eq. (5.1). 



33 



References 



[1] V. Silveira and A. Zee, Phys. Lett. B 161 (1985) 136; J. McDonald, Phys. Rev. D 50, 3637 
(1994); C. P. Burgess, M. Pospelov and T. ter Veldhuis, Nucl. Phys. B 619, 709 (2001); 
C. E. Yaguna, JCAP 0903 (2009) 003; M. Farina, D. Pappadopulo and A. Strumia, Phys. 
Lett. B 688 (2010) 329; W. -L. Guo, Y. -L. Wu, JHEP 1010, 083 (2010); Y. Mambrini, 
Phys. Rev. D 84 (2011) 115017; A. Djouadi, O. Lebedev, Y. Mambrini and J. Quevillon, 
Phys. Lett. B 709 (2012) 65. 

[2] K. Agashe, R. Contino and A. Pomarol, Nud. Phys. B 719, 165 (2005). 

[3] B. Gripaios, A. Pomarol, F. Riva, J. Serra, JHEP 0904 (2009) 070. 

[4] J. Galloway, J. A. Evans, M. A. Luty and R. A. Tacchi, JHEP 1010 (2010) 086. 

[5] T. A. Ryttov and F. Sannino, Phys. Rev. D 78 (2008) 115010. 

[6] Y. Hosotani, P. Ko and M. Tanaka, Phys. Lett. B 680 (2009) 179 . 

[7] J. L. Diaz-Cruz, Phys. Rev. Lett. 100 (2008) 221802. 

[8] G. F. Giudice, C. Grojean, A. Pomarol, R. Rattazzi, JHEP 0706 (2007) 045. 

[9] P. Gondolo and G. Gelmini, Nucl. Phys. B 360 (1991) 145. 
[10] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011). 
[11] C. E. Yaguna, JHEP 1108 (2011) 060. 

[12] M. Frigerio, T. Hambye and E. Masso, Phys. Rev. X 1 (2011) 021026. 

[13] LEP Higgs Working Group for Higgs boson searches and ALEPH and DELPHI and 
CERN-L3 and OPAL Collaborations, [hep-ex/0107032]. 

[14] ATLAS Collaboration, arXiv: 1202. 1408 [hep-ex]; S. Chatrchyan et al. [CMS Collabora- 
tion], arXiv: 1202. 1488 [hep-ex]. See also the recent slides at http://indico.in2p3.fr/ 
conf erenceDisplay . py?conf Id=6001. 

[15] A. Falkowski, Phys. Rev. D 77, 055018 (2008); I. Low and A. Vichi, Phys. Rev. D 84, 
045019 (2011); A. Azatov and J. Galloway, arXiv:1110.5646 [hep-ph]. 

[16] E. Aprile et al. [XENONIOO Collaboration], Phys. Rev. Lett. 107, 131302 (2011). 

[17] R. Bernabei, P. Belh, F. Cappella, R. Cerulh, C. J. Dai, A. d'Angelo, H. L. He and A. In- 
cicchitti et a/., Eur. Phys. J. C 67, 39 (2010); R. Bernabei et al. [DAMA Collaboration], 
Eur. Phys. J. C 56, 333 (2008). 



34 



[18] C. E. Aalseth, P. S. Barbeau, J. Colaresi, J. I. Collar, J. Diaz Leon, J. E. Fast, N. Fields 
and T. W. Hossbach et al, Phys. Rev. Lett. 107, 141301 (2011). C. E. Aalseth et al. 
[CoGeNT Collaboration], Phys. Rev. Lett. 106, 131301 (2011). 

[19] G. Angloher, M. Bauer, I. Bavykina, A. Bcnto, C. Bucci, C. Ciemniak, G. Deuter and 
F. von Feilitzsch et al, arXiv:1109.0702 [astro-ph.CO]. 

[20] D. Hooper, arXiv: 1201. 1303 [astro-ph.CO]; 

C. Kelso, D. Hooper and M. R. Buckley, arXiv: 1110.5338 [astro-ph.CO]. 

[21] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996). 

[22] M. Cirelh, G. CorceUa, A. Hektor, G. Hutsi, M. Kadastik, P. Panci, M. Raidal and 
F. Sala et al, JCAP 1103, 051 (2011); P. Ciafaloni, D. ComeUi, A. Riotto, F. Sala, 
A. Strumia and A. Urbano, JCAP 1103, 019 (2011). 

[23] J. Lavalle, Phys. Rev. D 82, 081302 (2010). 

[24] O. Adriani et al [PAMELA Collaboration], Phys. Rev. Lett. 105, 121101 (2010) . 

[25] S. Galh, F. locco, G. Bcrtone and A. Melchiorri, Phys. Rev. D 80 (2009) 023505 
[arXiv:0905.0003 [astro-ph.CO]]. 

[26] D. Feldman, Z. Liu and P. Nath, Phys. Rev. D 79 (2009) 063509. 

[27] M. Ibe, H. Murayama and T. T. Yanagida, Phys. Rev. D 79, 095009 (2009) . 

[28] W. -L. Guo and Y. -L. Wu, Phys. Rev. D 79, 055012 (2009) . 

[29] M. Ackermann et al [LAT Collaboration], Photon Spectrum," arXiv:1205.2739 [astro- 
ph.HE]. 

[30] P. J. Fox, R. Harnik, J. Kopp and Y. Tsai, arXiv: 1109.4398 [hep-ph]. 

[31] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, Phys. 
Rev. D 82 (2010) 116010; K. Cheung, P. -Y. Tseng, Y. -L. S. Tsai and T. -C. Yuan, 
arXiv:1201.3402 [hep-ph]. 

[32] J. R. Espinosa, C. Grojean, M. Muhlleitner, JHEP 1005 (2010) 065; D. Carmi, 

A. Falkowski, E. Kuflik and T. Volansky, arXiv:1202.3144 [hep-ph]; A. Azatov, R. Con- 
tino and J. Galloway, arXiv: 1202.3415 [hep-ph]; J. R. Espinosa, C. Grojean, M. Muh- 
lleitner and M. Trott, arXiv:1202.3697 [hep-ph]. 

[33] J. R. Espinosa, T. Konstandin and F. Riva, Nucl. Phys. B 854 (2012) 592 ; J. R. Espinosa, 

B. Gripaios, T. Konstandin and F. Riva, JCAP 1201 (2012) 012 . 

[34] J. Kehayias and S. Profumo, JCAP 1003 (2010) 003 . 



35 



[35] D. B. Kaplan, Nucl. Phys. B 365 (1991) 259. 

[36] C. G. Callan, Jr., S. R. Coleman, J. Wess and B. Zumino, Phys. Rev. 177 (1969) 2247. 

[37] J. Mrazek, A. Pomarol, R. Rattazzi, M. Redi, J. Serra and A. Wulzer, Nucl. Phys. B 
853 (2011) 1. 

[38] K. Agashc, R. Contino, L. Da Rold and A. Pomarol, Phys. Lett. B 641, 62 (2006) . 

[39] R. Contino, L. Da Rold and A. Pomarol, Phys. Rev. D 75 (2007) 055014 . 

[40] H. -Y. Cheng and C. -W. Chiang, arXiv: 1202. 1292 [hep-ph]. 



36 



